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We study a quantum mechanical system consisting of up to three identical dipoles confined to move 
along a helical shaped trap. The long-range interactions between particles confined to move in this 
one dimension leads to an interesting effective two-particle potential with an oscillating behaviour. 

For this system we calculate the spectrum and the wave functions of the bound states. The full 
quantum solutions show clear imprints of the tendency for the system to form chains of dipoles 
along the helix, i.e. a configuration in which the dipoles are sitting approximately one winding of 
the helix apart so that they can take maximal advantage of the strong head-to-tail attraction that 
is a generic feature of the dipole-dipole interaction. 


I. INTRODUCTION 

Cold atomic gas physics has advanced to the stage at 
which simulation and exploration of quantum phenom¬ 
ena can now be done at unprecedented levels of accuracy 
and with great potential for engineering systems that are 
otherwise very hard to access in other fields of research 
[IHi. An important recent advance is the ability to cre¬ 
ate homo- and heteronuclear molecules or trap atoms 
with large permanent magnetic moments at low temper¬ 
atures which has many interesting applications |6l-ll4|. 
The head-to-tail attraction between such molecules with 
dipole moments is a concern, as it limits the timescales 
at which an experiment can operate due to strong loses. 
One way to solve this is to confine the molecules in lower 
dimensional traps, where the head-to-tail attraction can 
be suppressed [lafig. 

The long-range dipole-dipole interaction has a very in¬ 
teresting interplay with the geometry of the lower dimen¬ 
sional systems. Here we consider the one dimensional 
geometry of a helix. This kind of system can be real¬ 
ized with Laguerre-Gaussian beams with nonzero angu¬ 
lar momentum [l7l - [^ or by trapping the atoms in the 
evanescent field surrounding an optical fiber The 

detailed interplay between the non-trivial long-range in¬ 
teraction of both charged ions and dipolar particles on a 
helical geometry has given rise to a number of theoretical 
studies into these setups [28l - [^ . 

In this work we consider dipolar particles on a single 
helix which confines the movement of the dipoles. We will 
be assuming that the helical trap is strongly confined in 
the directions that are locally perpendicular to the helix 
such that we may ignore the motion of the particles in 
these transverse directions. In this setup we investigate 
the formation of quantum mechanical bound states of two 
and three particles on the helix, we study their structure 
and the criteria for formation. The paper is organized 
as follows. In section II we introduce the helix and the 
dipole-dipole interaction. In section III we look at two 
dipoles on a helix, we calculate the bound states and 
look at their energy and their size. In section IV we 
introduce another dipole and look at the bound states of 
three dipoles on a helix. We conclude in section V. 


II. THE HELIX MODEL 

A helix is defined by two parameters, the radius R, and 
the pitch h. The radius defines the circumference while 
the pitch defines how far one moves along the symmetry 
axis of the helix during one revolution. The geometry 
of the helical setup is illutrated in figure [TJ To describe 
the position along the helix it is convenient to use the 
arc length s measured along the helix, but for numerical 
purposes we will instead use the angle around the central 

symmetry axis 4> = s/a, where a = ^R? + (^)^- It 
worth noting, that this angle (j) is not limited to be smaller 
than 27r, anything more than that just means more than 
one revolution on the helix. The relation between the 
coordinates on the helix, and the Cartesian coordinates 
in three dimensional space is 

{x,y,z) = {Rsm(l),Rcos(j),h-^), (1) 

Ztt 

where R is positive and </> runs through an interval be¬ 
tween (/min and (/max ■ On the helix we place between one 
and three dipoles, all of equal mass m, and dipole mo¬ 
ment d, all the dipoles are aligned by an external field. 

The long-range interaction of two particles on the helix 
are given by the usual dipole-dipole interaction between 
two dipoles in three-dimensional space. The two-body 
interaction potential, V, between two dipoles at position 
Vi and Tj is 

y{ri,rj) = [d- ti-3(d-r) (d- r)], (2) 

where r = — rj\ is the distance between the dipoles, 

and f = {ri — rj)/r is the unit vector in the direction 
connecting the two dipoles. We use the unit where 
eo is the permitivity, corresponding to electric dipoles, 
but in principle the formalism is applicable for magnetic 
dipoles if we use the unit where /xg is the vacuum 
permeability. Here we will assume that an external field 
is used to align the dipole moments of all particles along 
the symmetry axis which we take to be the z-axis (see 
figure [T]). As mentioned above, instead of describing the 
position of the dipoles by their Cartesian coordinates it 
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FIG. 1. A representation of a helix of radius R and pitch 
h, with dipoles positioned at si and S 2 , respectively. Here s 
is the arc length measured along the helix. We assume that 
an external field is applied which directs the dipole moments 
along the z-axis. The cartesian axes are indicated on the left 
at the bottom. 


turns out to be easier to use their position along the helix. 
This can be done through the coordinate transformation 
in Eq. ([T]). If the position of the two dipoles along the 
helix are <j)i and tpj the two-dipole potential in Eq. (ED 
becomes 




dtreo 


( 3 ) 


2R^ [1 - cos {(j>i - (j)j)] - 2h‘^ - (j)j)/ (27r))^ 

^2i?2 [1 - cos - (j)j)/ (27r))^^ 


We see that for dipoles trapped on the helix with dipole 
moments pointing in the the ^-direction the two-dipole 
potential in Eq. ED only depends on the distance between 
the two dipoles (f> = (t>i — (l)j measured along the helix. 

The two-dipole potential as a function of this dipole 
separation, is shown in figure [21 The long range part 
of the potential falls of as like the distance depen¬ 
dence from Eq. Q in three dimensions. But the poten¬ 
tial also contains oscillations at smaller distances. They 
disappear when ^ > 1 + (27r)^-^. The corresponding 
minima in the potential arise whenever the dipoles are 
above each other. It is worth noting that the minima are 
not exactly at multipla of 27r, which would correspond to 
the dipoles being directly above each other. Instead the 
minima are located at a distance from each other that is 
slightly smaller than 27r as it is a compromise between 
optimizing the angle, and minimizing the three dimen¬ 
sional distance between the dipoles. 

When two dipolar particles come very close to each 
other we expect that strong interactions will cause for 
instance chemical reactions and losses from our setup. 
This is not the regime we are interested in and we there¬ 
fore look for a parameter regime in which the short-range 
behavior of the dipole-dipole potential is purely repulsive. 



FIG. 2. The reduced potential, V = A.-k eoV / cR , of two 
dipoles as a function of the relative angle, (f> = <f)i — 4'] t ssP" 
arating the particles on the helix. The helix parameters are 
chosen to be = it (blue solid) and h = 1.6it (black dashed). 


In this regime, we may safely ignore any chemical reac¬ 
tions and strong losses. This requirement puts limita¬ 
tions on the geometry of the helix. If -^ > v^tt, then the 
two-dipole potential becomes attractive, and would allow 
distances where f.x. van der Waals forces would domi¬ 
nate. We thus restrict ourselves to cases where ^ < v^tt. 


III. THE TWO-BODY PROBLEM 


We now place two dipoles on the helix. We assume 
that they are only subject to their mutual interaction 
which only depend on the relative distance between the 
dipoles as in Eq.(I3D. It is convenient to make a change 
of variables, from the positions of the two dipoles (j)i and 
(j) 2 , to relative angular distance (j) = (j)i ~ ^ 2 , and center 
of mass ^ = (j>i + (j> 2 - A straightforward if somewhat te¬ 
dious calculation now shows that the relative and center 
of mass motion decouples in the Hamiltonian. The center 
of mass motion will be governed by a free particle Hamil¬ 
tonian which can be trivially solved for given boundary 
conditions. In this paper we will ignore the center of 
mass motion which is not important for the question of 
bound state formation. The Hamiltonian governing the 
relative motion of the two dipoles takes the form 


H = 


2)10^ d(jP' 


+ V{<l2) 


( 4 ) 


Here ^ = m/2 is the reduced mass of the two dipoles. 
Note that one must be careful with quantization in ge¬ 
ometries with non-trivial curvature. However, the regular 
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helix we consider in this paper has the property that its 
curvature is in fact constant and one may thus do the 
transformation from Cartesian to curvilinear coordinates 
in the standard manner [s^, . It is convenient to write 

the Hamiltonian in natural units, iJ, which has the form 


where /3 is the potential strength given by 

/ rv \ 2 


(!)■ 


2t:€oRTP‘ 

and V is the reduced potential where 
l-cos<^- 


V(c/>) = 


2 [1 - cos(j)] + 


5/2 ■ 


(5) 


( 6 ) 


(7) 


We thus measure lengths in units of a and energies in 
units of Because —>• oo as 0 —>■ 0, the 

wavefunction tplcj)) has to be zero at </) = 0. This means, 
that we only need to calculate the wavefunction for </> > 0, 
and then construct the wavefunction to be either odd or 
even depending on whether dipoles are identical fermions 
or bosons. 


A. Solutions to the two-dipole system 

We solve the Schrddinger equation to obtain the wave 
function, tp, for the relative motion of the two dipoles in 
a box of size (p G [0 : 100], with closed boundary condi¬ 
tions, that is t^(0) = ■!/)(100) = 0. Since the dipole-dipole 
interaction is repulsive at short range the wave function 
must go to zero as the two particles coincide in space. We 
may therefore reduce the complexity of the problem and 
solve in the region (pi > p 2 only (or vice versa). The full 
wave function can subsequently be found by extended to 
the opposite domain using continuity and considering the 
statistics of the particles (even function for bosons and 
uneven function for fermions). We start out by solving 
it for the case of /3 = 1 and -^ = 1. In figure O the 
wave functions of the four lowest eigenstates are shown. 
They all have a negative energy corresponding to bound 
states of two dipoles on a helix. The solid red curve is 
the ground state. It has a single maximum, aX (p k. 2'k 
corresponding to where the minimum in the two-dipole 
potential is located. The first excited state is the dashed 
green curve, it has a single node close to p = 10, between 
the first two minima of the two-dipole potential. The 
maximum is close to the second minimum in the poten¬ 
tial corresponding to the two dipoles being separated by 
two windings on the helix. The short-dashed blue curve 
is the wavefunction of the second excited state. It has 
two nodes between the first three minima of the poten¬ 
tial, and corresponding maxima at the first two minima. 



FIG. 3. The wave functions for four lowest energies with /3 = 
1. Note how the ground state ipo has its peak at p = , 

the classical equilibrium distance. 


but because it has to be orthogonal to the two lower ly¬ 
ing states it has a large part beyond the third minimum. 
The last curve in dotted purple is the wavefunction of 
the third excited state. As seen, the curvature at large p 
is negative, whereas the other three had a positive cur¬ 
vature. This is because of the boundary condition at 
p = 100 where the wavefunction is forced to be zero for 
our calculations. We will return to this boundary effect 
below as it will also show up in our calculations of the 
case with three dipoles on the helix. 

In the case of /3 = 1 shown in figure O there are three 
bound states. In figure 01 the energies of the four lowest 
states of two dipoles are shown. In the figure one can 
see that the number of bound states increases with the 
potential strength as expected. Whether a bound state 
exists for all positive values of (3 is not possible to answer 
a priori. Generally in one dimension, if the integral over 
the potential is negative, that is if / V{p)dp < 0 then a 
bound state exists [3§ - l40l |. But in our case the poten¬ 
tial diverges at short distances and to calculate this inte¬ 
gral one would need knowledge of the before mentioned 
short-range forces that are out of the scope of this paper. 
The exact potential strength for which individual states 
become unbound is dependent on the numerical parame¬ 
ters, but the scaling on the number of bound states with 
/3 is always seen. 


B. The Size of the bound states 

In this section we explore how a ground state solution 
behaves, when the dipole moment changes. Changing the 
dipole moment corresponds to changing the dimension¬ 
less parameter /?. For a large dipole strength (/3 large), 
we expect the system to resemble a harmonic oscillator. 
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FIG. 4. The lowest four eigenvalues of two dipoles in a helical 
trap for different potential strengths /3. The energy is in units 



with the two dipoles oscillating about the distance where 
the potential is at a minimum. The size of a harmonic 
oscillator follows the simple relation 

< (j)^ >(x - (X ^ (8) 


In this case the harmonic oscillator is not centered around 
0 but instead shifted by the equilibrium distance (/)o so 
accordingly the value of < > should be shifted by ^q, 

where (j)o is where ^ = 0. So we expect < > to scale 

as 


< (j)^ >= Cl — +C2(j}, 


(9) 


where ci and C 2 are fitting parameters. We are not inter¬ 
ested in their specific value. In figure [S] the size < cj)^ > 
of the ground state is shown as a function of the param¬ 
eter P, while the solid red line is a fit to Eq. (jQ)). One 
sees that in the large P limit the size of the ground state 
behaves as that of a harmonic oscillator. 

For small P and thus small energies, we expect the sys¬ 
tem to approach free particles, in the sense that we expect 
the wave functions to asymptotically approach the form 
= e~'^^ where k = y/^. Using the asymptotic 
form of the wave function we can calculate the size of the 
bound states as a function of either the energy E, or the 
strength of the potential p. Because the wave function 
for the asymptotic behaviour is not normalized we can 
calculate {(fP) using 




(jPe 

/o°° 


1 


( 10 ) 


So for low energies we expect {p"^) to scale as -g according 
to equation (fTUl) . In figure [5] we show the product of size 



P 

FIG. 5. The expectation value of p^ as a function of /3 is 
shown with plus symbols. The solid line displays a fit of the 
data to the functional form given in Eq. 




FIG. 6. The expectation value of p^ as a function of binding 
energy, for the two-body ground state. 


of the ground state {p"^) — pQ and the ground state energy 
E as a function of the energy. 

If the scaling follows the we expect this to approach 
a constant as E decreases. In the figure there is a region 
where it is constant around E = —0.015. However, for 
very small energies the product increases rapidly. This is 
because of boundary effects where the size of the ground 
state approaches the size of the box in which we solve the 
problem numerically. Alternatively, we see that for the 
helix sizes we work with the ground state for very small 
energies becomes so large that the asymptotic behavior 
is not reached. However, the state does push the bound¬ 
aries of the geometry as it attempts to extend toward 
large distances. We have run numerical tests to confirm 
that once the wave function becomes very broad at low 
energy, its size scales with the size of the bounding box. 
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IV. THE THREE-BODY PROBLEM 


We now add a third dipole to the helix. It is of the 
same type as the other two, that is the same mass and 
dipole moment. They all interact through the two-dipole 
interaction in Eq. ©• The position of each of the dipoles 
is denoted 4’i^4’2 and ^ 3 . As is typically the case, it is of 
course very convenient to separate out the center of mass 
motion of the three dipoles. However, one must keep in 
mind that we are considering the non-trivial geometry of 
the helix and it requires some consideration to work this 
out. The criterion is that all the particles have the same 
mass, m, which is the case we are interested in here. The 
separation can now be done as in other equal-mass three- 
body systems through the coordinate transformation 


a; = ^ ((/>! - (/> 2 ) 

2 = -^ (</>!+ </>2 + h ) ■ 


( 11 ) 


With no external potential along the helix, the three 
body problem now reduces to a two dimensional prob¬ 
lem, in the two relative coordinates x and y as the z- 
coordinate only describes a center of mass motion. Note 
that the three coordinates x, y, and 2 : are atypical in the 
sense that they are built from angular variables and thus 
are intrinsically dimensionless. Without an external po¬ 
tential along the helix, this center of mass motion can be 
easily worked out. We use open boundary conditions on 
the helix and thus require the wave functions to vanish 
on the top and bottom of the helix which implies that 
the center of mass behaves simply as a particle in a box. 

The two dimensional Schrodinger equation that de¬ 
scribes the relative motion of the three dipoles should 
contain five terms, two terms from the kinetic energy of 
the relative motion between the dipoles, and three po¬ 
tential terms from the pairwise dipole-dipole interaction 

Here V is the two-dipole potential of Eq.d?]). All terms 
are written in natural units using /3 as in the case of 
two dipoles above. The Hamiltonian is also in units of 
yc? as in Eq. ©■ Because the two-dipole potential 
is repulsive at short range, we use the same trick as for 
two dipoles to reduce the complexity of the system, we 
only solve it for ipi > (j)2 > 'Ps- then depending 

on whether the dipoles are bosons or fermions, the full 
symmetric or antisymmetric wave function can be recon¬ 
structed for all values of 4>i,4>2 and <^ 3 . The restriction 
to 4>i > 4>2 > 4>3 corresponds to a restriction on x and y, 


'ip{x,y) 
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FIG. 7. The ground state wave function of three dipoles on 
a helix with a potential strength /3 = 1. The dominant peak 
corresponds to a configuration where the dipoles are sitting 
(approximately) one winding apart. Notice that the x and y 
coordinates are combinations of angles and thus are dimen¬ 
sionless (see Eq. (HU). 


namely that y > -^x. For definiteness, we will be work¬ 
ing with a dipolar strength parameter of ,3 = 1 when 
we plot the wave functions (below we will explore varia¬ 
tions with /3 and its effect on the relative distances within 
the ground state). In addition, we take h = R, i.e. the 
pitch and the radius are the same. The two-dimensional 
wave equation is solved on a grid with a; € [0 : 160] 
and y > We observe the same behavior as for two 

dipoles with regards to the size of the states which also in¬ 
creases with decreasing /?. From a practical point of view 
this means that for smaller /3 the size of the grid has to 
be increased which slow down the numerical solution and 
reduces accuracy. We therefore restrict ourselves to the 
case of /? = 1 for three dipoles. 


A. Properties of the solutions 

The ground state wave function of the three-body 
problem is shown in figure [3 in a three-dimensional level 
plot that includes a two-dimensional projection onto a 
contour plot for clarity. The ground state displays a 
prominent single peak at around x ~ 4.4 and y ~ 7.7. 
This peak corresponds to a configuration where the three 
dipoles are placed directly above each other only sepa¬ 
rated by (approximately) one winding. More precisely, 
it corresponds to a configuration where (<^ 1 , (() 2 , (/is) ^ 
( 47 r, 27r, 0). We thus see a very nice consistency with re¬ 
spect to the case of two dipoles in which we also find 
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a dominant configuration with the dipoles placed about 
one winding of the helix apart in order to exploit the 
attractive interaction in the head-to-tail setup. 

We now consider the low-energy spectrum of three 
dipoles on a helix by representing the wave functions in 
contour plots as functions of x and y coordinates. This is 
shown in figure HI As discussed above, the wave equation 
has been solved only for i/Sj/ > x. Solid (black) lines in¬ 
dicate where x = Ta/Sj/ in all the panels in figure HI To 
extend the results to the full (x, j/)-coordinate space, one 
needs to reflect the wave functions across the solid lines 
in the panels. Here we are assuming that the dipoles 
are identical bosons. This implies that the wave func¬ 
tion must be symmetric in the exchange of any pair of 
particles, and thus we must reflect across the solid lines 
and keep a positive sign. In the case of fermionic dipoles, 
the solution would be the same except for a sign change 
across the diagonal. 

The panels in figure HI show the first four states in the 
spectrum. In the upper left corner, we have the ground 
state which is identical to the state shown in a different 
manner in figure [T] It is a state with three head-to-tail 
dipoles, here extended across the solid lines so it is re¬ 
ally three copies of the central region (within the wedge 
traced by the solid lines containing the y = 0 axis). The 
first excited state in the top right corner has a much more 
detailed structure. It still contains a trace of the head-to- 
tail on successive windings structure of the ground state 
but now a lot of amplitude is move to larger values of 
X in the central region around y = 0 (within the cen¬ 
tral wedge). This corresponds to configurations where 
the one of the outer dipoles from the head-to-tail config¬ 
urations is now pushed one winding away from the two 
others (this can be done in two ways so we have y —^ —y 
symmetry here). Notice also the change in sign between 
the two regions with non-zero amplitudes. This is of 
course a result of the fact that higher excited states have 
additional nodes in the wave functions. We may relate 
this very directly to the physics we saw in the case of two 
dipoles in figure H] There we see that the first excited 
two-body state (long dashed (green) line) has a small bit 
of amplitude around angles of ^ ~ 27r but that most of 
its weight is around (f) ^ An and thus the two-body state 
has the dipoles sitting about two windings apart mainly. 
This is clearly also reflected in the contour plot of the 
first excited state for three dipoles seen in the top right 
panel of figure HI 

The second excited state of three dipoles is shown in 
the lower left panel in figure HI and has a structure that 
can now be simply understood given the two lower states. 
Its amplitude is dominated by pairs of dipoles which are 
now about two winding apart. In comparison to the first 
excited state it has zero amplitude of the simple head- 
to-tail configuration seen in the ground state. It is a bit 
harder to compare this to the two-body case as the sec¬ 
ond excited state for two dipoles has them sitting really 
far apart (at the level of (f> ^ 5n, see hgureHlli while the 
second excited state for three dipoles resembles more the 


physics in the first excited state in the two-body case only 
without the configurations where two of the dipoles are 
sitting one winding apart. The third excited state seen 
in the lower right panel of figure HI tells a similar story ex¬ 
cept that now the dipoles move even further apart as the 
amplitude is seen to move to larger values of x and y. We 
also notice that a small bit of amplitude comes back to 
the ground state configuration in the region near the ori¬ 
gin {x, y) = (0, 0). However, the third excited state is not 
fully converged as one can see by the lack of symmetry 
(up to a sign) for y ^ —y. It is the same sort of boundary 
effect that can be seen in figure H] for the third excited 
state of the two-body system. In spite of this numerical 
issue, the tendency of higher excited states should now 
be clear. In higher excited states the dipoles are pushed 
further and further away from each other and is consis¬ 
tent with the picture that we have from the two-body 
case. A nice feature is that the ground state configura¬ 
tion with three dipoles head-to-tail on three successive 
windings does indeed seems to make an appearance in 
higher excited states also so we do see that the first min¬ 
imum in the two-body potential in figure HI plays a very 
dominant role in this geometry. 

In order to further elucidate the configuration of the 
three dipoles of the helix we can consider the relative 
distances in angle between each pair of dipoles within the 
ground state. These distances are defined as (t)ij = (j)i — (j)j 
and we take the expectation value of this operator in 
the ground state. The results are shown in table U for 
different values of (3 with h = R. The distances are 
calculated in units of 27r. Because of the chosen ordering 
they are all positive, and ^13 = ^12 + ^!' 23 - As seen in the 
table, for larger values of /3 the ground state has a clear 
interpretation as a ’chain’ of three dipoles separated by 
one winding. This is seen to set in already for (3 = 1 and 
is accurate at the level of two decimal places already for 
[3 = 2. For smaller values of [3 we expect that kinetic 
terms will be more important and the dipoles would like 
to delocalize. In the table this is seen for a value of 
/3 = 0.25 where the expectation values of the distances 
are no longer close to 27r, i.e. the dipoles tend to be 
sitting further than one winding apart on average. In 
this regime of weak dipolar interactions the dipoles will 
tend to spread out to minimize kinetic energy while at the 
same time being able to take advantage of the attraction 
from several of the pockets seen in figure HI 


TABLE I. Relative angular distances in the ground state for 
three dipoles with h = R in units of 27r for different values of 

/3. 


p 

{<t>i2)l2n 

{(j>23)l2n 

{(j)i3)/2n 

P = 0.25 

1.50 

1.44 

2.94 

/3 = 1 

1.01 

1.01 

2.03 

P = 2 

1.00 

1.00 

2.00 
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FIG. 8. The wave functions of the four lowest states of three dipoles on a helix shown as contours. The bright (red) regions 
are positive while the dark (blue) region are negative. Notice that the x axis is vertical, while the y axis is horizontal. The 
solid (black) lines in panels indicate the lines where x = ±\/%y. These are the reflection lines of particle exchange discussed in 
the text. Here the particles are assumed to obey bosonic statistics (see the text). The ground state (top left), the first excited 
state (top right), the second excited state (bottom left), and the third excited state (bottom right) are shown. Note again that 
X and y are dimensionless as in figure [T] 


V. DISCUSSION AND OUTLOOK 

In the present paper we have considered the physics 
of dipolar particles that are confined to move on a one¬ 
dimensional helix. We first look at the dipole-dipole in¬ 
teraction on a helix, and how the interplay between the 
long range interaction and the peculiar geometry of the 
helix leads to a two particle potential with several min¬ 
ima of decreasing depth. These minima correspond to 
the attractive head-to-tail configuration and the decreas¬ 
ing depth of successive minima is a result of the dipoles 
being an increasing number of winding of the helix apart. 
Our main question concerns the formation of bound state 
in this non-trivial system. In particular, what the forma¬ 
tion criteria for two-body bound states are and what type 
of bound states form with more dipoles. As our main fo¬ 
cus we use the case of three dipoles. 

In the strong interaction limit both two- and three- 
body bound states correspond to the dipoles being an 
integer number of windings apart in the quantum ground 
state of the system. However, for weaker interactions the 
bound states increase in size, and pairs of dipoles can 
no longer be said to be a certain number of windings 
apart as they become effectively delocalized across sev¬ 
eral windings of the helix. For three dipoles we show 
that for moderate dipole strengths they form a short yet 
well-defined chain of three dipoles sitting immediately 


underneath each other. 

This affinity for chain formation was discussed previ¬ 
ously in the limit of very strong interactions where classi¬ 
cal crystal formation on the helix is expected . There 
it was shown how the dipoles would form chains up the 
helix where each dipole was approximately one winding 
away from its neighbors on either side. In one dimen¬ 
sion such long range orders are not possible in the quan¬ 
tum regime, and one would instead expect the formation 
of a Luttinger Imuid in such a system as has been dis¬ 
cussed in Ref. |28j . This should be more pronounced 
for smaller values of j3 where kinetic terms are sizable 
and the particles will tend to delocalize, i.e. not merely 
stay approximately fixed in the minima provided by the 
attractive head-to-tail configuration. The formation of 
chains in related geometries in both two- [dll, and 
one-dimensional (^.1^ setups, and this chain formation 
is expected to persist and be in the many-body case also 
in the quantum regime [dl|, . The strongly interact¬ 

ing regime can be explored by using harmonic approxi¬ 
mations [d^ to the full dipolar interaction for instance 
in the study of the thermodynamic properties of dipolar 
chains [d^- As we discussed above, a harmonic behavior 
is also seen in the helical geometry for large /3 which is 
merely a reflection of the fact that the dipolar potential 
allows such an approximation for large dipole moments 
in any geometry where the head-to-tail configuration is 
possible. In future studies it would be interesting to ex- 












tend the system to slightly longer chains and study the 
thermodynamics for instance using the harmonic approx¬ 
imation for strong interactions. 
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